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A LOCAL ANALYSIS OF THE AXI-SYMMETRIC NAVIER-STOKES FLOW 
NEAR A SADDLE POINT AND NO-SLIP FLAT BOUNDARY 

PEN-YUAN HSU, HIROFUMI NOTSU, AND TSUYOSHI YONEDA 


Abstract. As one of the violent flow, tornadoes occur in many place of the world. In order 
to reduce human losses and material damage caused by tornadoes, there are many research 
methods. O ne of the effect ive methods is numerical simulations such as the work in a recent 
article Ushihar a et al\ (|20lD . The swirling structure is significant both in mathematical analysis 
and the numerical simulations of tornado. In this paper, we try to clarify the swirling structure. 
More precisely, we do numerical computations on axi-symmetric Navier-Stokes flows with no¬ 
slip flat boundary. We compare a hyperbolic flow with swirl and one without swirl and observe 
that the following phenomenons occur only in the swirl case: The distance between the point 
providing the maximum velocity magnitude |t| and the z -axis is drastically changing around 
some time (which we call it turning point). An “increasing velocity phenomenon” occurs near 
the boundary and the maximum value of |t| is obtained near the axis of symmetry and the 
boundary when time is close to the turning point. 


1. Introduction 


In this paper we consider a local behavior of the 3D-Navier-Stokes flow near a saddle point 
(with hyperbolic flow configuration) and no-slip flat boundary. The Navier-Stokes equations 
with no-slip flat boundary are expressed as 

(1.1) dtv + (v • V)v — isAv + Vp = 0 in x [0,T), 

v 0 = vjt= 0 , ?Tm 3_ = o, V-v = 0 in M+x[0,T), 

where v is a vector field representing velocity of the fluid, and p is the pressure. The term 
“hyperbolic flow configuration” which used here and after means there is 5 > 0 (depending on t) 
such that v(t, x) • e z > 0, v(t, x) • e r (x) < 0, or v(t, x) • e z < 0, v(t, x) • e r (x) > 0 for 0 < \xh\ < 5 
and 0 < £3 < 5, where e z — (0,0,1), e r (x) = {x\/\xh\,X 2 /\xh\,fy and \x^\ = \Jx\ + x\. For the 
definition of 2D-hyperbolic flow configuration, we just need a simple modification. 


Recently, the Euler flow with a sac 


die point (with hyperbolic flow configuration) has been 


extensively studied. iBourgain fc Lil (2 013 ) obtained strong local ill-posedness results in the 


Sobolev spaces W n / p+1,p for any 1 < p < 00 and in the Besov spaces Bp(q +1 with 1 < p < 00 
and 1 < q < 00 and n — 2 or 3 by using a combination of Lagrangian and Eulerian techniques 
with a saddle point (with hyperbolic flow confi guration) structure. I n par ticular , they settled 
the borderline Sobolev case H n / 2+1 . After that lElgindi fc Masmoudil (120141 ) and Bourgain fe Li 


(2015) producec 
( 2014 ) (see also 


similar re s ults i n C 1 case (also C m case). On the ot 


Itoh et al 


( 20141 ): Kiselev & Zlatos (2Q14j); 


Xu 


EoiT 


Kiselev & Sverak 


rer hand, 


for related topics) showed 


a two-dimensional Euler flow with a saddle point (with hyperbolic flow configuration) in a disk 
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for which the gradient of vorticity exhibits double exponential growth in time for all times. Their 
estimate is known to be sharp, namely, the double exponential growth is the fastest possible 
growth rate. These results show that the saddle point (with hyperbolic flow configuration) bring 
(some kind of) unstabilizing effects. 

Now let us look back the history of Navier-Stokes equations briefly. Modern regularity t heory 
for the solutions to the Navier-Stokes equations began with the works of lLeravl (11934) and IHopf 
(1951). They showed the existence of a weak solution v : [0, oo) x K 3 A M 3 lying in the class of 
L°°(0, oo; L 2 (M 3 )) n L 2 (0, oo; ij 1 (M 3 )) which satisfies the global energy inequality, where H 1 is 
a homogeneo us Sobolev sp a ce with degre e one. _ 


After that, 


Prodi 


) lev sp ace witri 
(1959), Serrinl 


(119631 ). iLadvzhenskaval (119671 ). and their joint efforts lead to 


the following Prodi-Serrin-Ladyzhenskaya criterion for the Leray-Hopf weak solutions. 

Theorem 1.1 (Prodi-Serrin-Ladyzhenskaya). Let v G L°°(0, oo; L 2 (R 3 )) H L 2 (0, oo; iL 1 (M 3 )) be 
a Leray-Hopf weak solution to (EUX which also satisfies v G L p (0, T; L q (R 3 )), for some p , q 
satisfying | + | = 1, with q > 3. Then, the solution v is smooth and unique on (0,T] x R 3 . 

It is also worthwhile to men tion th at t h e exc eptional case of v G L°°(L 3 ) was finally es¬ 
tablished in the work bv llEscauriaza et a/J (2003). After the appearance of the Prodi-Serrin- 
Ladyzhenskaya criterion, many different regularity criteria of solutions to HD was established 


b r ex ample, a regularity criterion 


Chan & Yoneda ( 201 2)). Besides 


by many researchers working in the regularity theory of (HH). 
along streamlines (characteristic curves) was con structed (see 
these, other important works such as lGiga et all (2014), in which type I blow up was excluded 
for solutions to (11.11) under a regularity condition on th e vo rticity direction in the half space to 
the case of the no-slip boundary conditon 
pioneer work in this field), and 


Chen et al 


see also 


(2009); 


Constantin <fe Feffermanl (119931 ) , which is the 


Koch et al. 


(I 2 OO 9 I ) for axisymmetric solutions 


to (jl.ip (the equation to the axi-symmetric case can be rewritten to (jl.3j) ) in the whole space, 
are attracting a lot of attentions. We now recall that a singularity of a Navier-Stokes solution v 
at a time T is called Type I if 

C 


( 1 . 2 ) 


sup \v(x, t) | < 


Vr^t 

for some C > 0. IfMhe_singularity of v does not satisfy the condition HI, we say Type II 


Koch et al 


singularity (see 


equations is expressed as 


( 2009! ) for the type of singularity). The axi-symmetric Navier-Stokes 


u 


(1.3) 


d-fcUf —(— Uj'dj'Uj' u z d z u r —|— d r p — 

r 

0,0 | r\ . U r UQ _ 

OtUQ T u r o r UQ T u z o z uq T- — 

r 

dtu z ~\~ u r d r u z + u z d z u z d z p = 
d r (ru r ) 

-b d z u z = 


A Ur 

’ 

A U ° 

^UQ 2 5 


A u z , 

0, 
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where u r = u r (r,z,t), uq = ue(r,z,t) u z = u z (r,z,t) and A = d% + (l/r)<9 r + d z . The vector 


valued function v := u r e r + uqcq + with e r := 


Oi,x 2 ,0), e$ := 


and = (0, 0,1) represents velocity of the fluid, and p is the pressure. 




I t is known that ax i -symmetric solutions with no swirl (namely, up = 0 ) have to be regular. 


See 


Koch et al 


Caffarelli et al 


Ladvzhenskaval (119681 ) ; lUkhovskii fc Iudovichl ([1968) . On the other hand, 


(1982) showed that the axi-symmetric Navier-Stokes solution could only blow 


up on the axis. 
In Choi et al. 


( 20141 ). they considered a ID transport equation with an additional variable 
which comes from a square of the swirl component uq of the velocity field in the 3D axi-symmetric 
Euler flow. They showed that the transport equation can exhibit finite-time blow-up from 


smooth initial data using a contradiction argument, and their conclusion in Choi et al. (20141 ) 
is that in the 3D axi-symmetric Euler flow with the flat boundary, the best chance for possible 
singularity formation seems t o b e at the saddle point (at the axis of symm etry intersect the 
boundary) (see also Luo & Hou (2014))). On the other hand, in IKand (2004), he constructed a 
regularity theory. More precisely, suitable weak solutions of the three-dimensional axisymmetric 
Navier-Stokes equaions in a half space with no swirl case are Holder continuous up to the 
boundary exept for the origin. His result suggests that even in the Navier-Stokes flow case, the 
best chance for possible singularity formation may be at the saddle point. 

Although there are many fruitful results based on mathematical analysis as we recalled above, 
it is not easy to analyze locally such fluid mechanics to go a step further mathematically. Thus, 
it should be effective for us to attempt numerical approach. 

Besides, from the above mathematical literature, the saddle point seems to be a key place 
and the hyperbolic flow with swirl might be a ke y structure. Actua ll y, the swirl ra t ios ar e 
very s i gnifi cant in t he researc hes o 


tornadoes, see IWan fe Ding! (2005); 


Ishihara et al. 


( 2011 ); 


Nolan (2012); Ishihara & Liu (2014). There are several methods for researching tornadoes. Al¬ 


though the studies on real tornadoes and simulated tornadoes in laboratory are also important, 
numerical simulations provide a safe and cost effective way to analyze the behavier of torna¬ 
does. We concentrate on the numerical computation in the swirl case and the behavior near 
the axis of symmetry (saddle point) at the boundary. From a mathematical point of view, the 
saddle point bring some kind of unstabilzing effects. In the point of view from the studies of 
tornadoes, the behavior insides the core of tornado (the center region near the z-axis) is signif¬ 
icant and might be very different to it outsides the core. The behavior near the ground (lower 
boundary) is also very significant for preventing the damage caused by tornadoes. Our numer- 
ical result can be compa r ed wi t h the mathematica l literature (especially regularity results, see 


Constantin fc Feffermanl (119931 ): iGiga et all (120141 )) and also with the two-celled vortex struc¬ 


ture in the tornado researches (especially numerical simulations). For more references, refer to 
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Giga et all (|2014l ) for regularity results, refer to lLuo fc Houl (120141) for numerica l studi es of the 


Navier-Stokes and Euler equations, and refer to iNolanT ( 20121 ); 
studies of tornado-like vortices. 

In the next section, we will show the numerical computation on the difference between the 
swirl case and the no swirl case. We also show that the maximum value of \v\ occurs at the 
place near the axis of symmetry (saddle point) and the boundary when the time approaches to 
the critical turning point (t=0.35 and t=1.0). (See figure 1,2,...,7, especially figure 7) 


Ishihara & Liu (120141 ) for the 


2. Setting of the initial data and numerical results for the axi-symmetric 

Navier-Stokes flow 

In this section we set an initial data of the axi-symmetric Navier-Stokes flow with a saddle 
point (with hyperbolic flow configuration). We will compare two flows: with swirl (uq ^ 0) and 
no swirl (uq = 0) cases. In our numerical computation, we use the following cylindrical domain 
Li: 

(2.1) Li := = (xi, £ 2 , z) G M 3 : —a < z < 4a, \jx\ + x\ < 1 

and impose no-slip boundary condition: 

(2.2) v — 0 on dLi. 

We set the initial data in the following manner. (See also figure 1 and remark 12.21 ) Let 
g?(a, e, cr) = (a 2 + e) a , we set the initial velocity for the swirl case as follows: 

(2.3) u z = cp(r,ei,-Pi)ip(z,e 2 ,-02), 

(2.4) p = g?(r, e 3 , —/3 3 )(^(z, e 4 , /? 4 ), 

(2.5) u r = sign (z)pu z , 

(2.6) uq = g>(r, e 5 , —/3 5 )g?(z, e 6 , —/3 6 ) 5 

where c and /% (i = 1,2,..., 6) are constants. In the following numerical calculation, we set all 
6i and ^ equal to 1. As for the no swirl case, we only change uq to zero. 

Remark 2.1. By our setting of the initial data, the initial velocity magnitude \v\ at (xi,x 2 ,z) = 
(0, 0, 0) (we call this point the center of the initial velocity) is larger than it in other places. 
Phenomenons observed in this paper are more clear when the center of the initial velocity is 
close to the lower boundary (but not on the boundary). This is the reason for our choice of 
computational domain Li like (12.1 j) instead of a symmetric one. See also appendix for comparison 
between the behaviors of two different centers. 


Let r > 0 be a time increment of the computation and h > 0 a (representative) mesh size. 
We perform the computation by a stabilized Lagrange-Galerkin (finite element) scheme ( iNotsul 
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2008 


Notsu fe Tabatal . 2008, 


2015a, 


b ) based on the idea 


(2.7a) 1 (x) - v k ~ l (x - v k ~ l (x)r)} - 2_ Av k + Vp k = 0, 

(2.7b) V • v k - 5 0 h 2 Ap k = 0, 


for k = 1, 2, • • • , where v(x, kr) is simply denoted by v k {x), Re is the Reynolds number, and the 
stabilization parameter is set as 5o = 1. We note that under some conditions the scheme has 
mathematical convergence properties of order 0{r + h) for the velocity in H l (yi) and for the 
pressure in L 2 (fl) and of order 0{r-\-h?) for the velocity in L 2 (fl). The maximum, minimum and 
average mesh sizes are 1.88 x 10 -2 , 1.48 x 10 -3 and 8.95 x 10 -3 , respectively, where the mesh size 
around z -axis is smaller than that of other part. The time increment is set as r — 1.25 x 10 -2 . 
In the following we show only numerical results by using the mesh size and the time increment 
mentioned above, where numerical results with a coarser mesh size and a larger time increment 
are qualitatively similar. See appendix for the discussion of the dependency of the numerical 
results on the discretization parameters h and t. 


Remark 2.2. The initial velocities for the swirl and no swirl cases do not satisfy divergence free 
condition and the no-slip boundary condition, the computational velocities after the first time 
step, however, satisfy both conditions numerically, where the former condition is satisfied in the 
sense that the equation (I2.7bl) holds. The structure of initial data is useful for observing the 
swirling flow. It is also reasonable for us to consider the initial data (I2.3I) - (I2.6I) as a pre-stage, the 
velocities after the first step as a real initial velocity in numerical computation. Although this 
construction of initial data is useful for observing phenomenons in numerical approach before 
the construction of more smooth initail data, a more careful choice of the initial data is desirable 
for mathematical analysis. 

Remark 2.3. The axial symmetry is not explicitly imposed in the three-dimensional computation 
by the stabilized Lagrange-Galerkin scheme, while the problem setting including the initial 
velocity has the symmetry. We perform the computation in order to show qualitative properties 
of the effect of the swirl. 


The following figures are our numerical results for a = 1/8 and Re = 1,000, 5,000, 10,000 
and 50,000. Hereafter, we call the point attaining the maximum value of |t| “the maximum 
point of \v\” and the cross-section {x E Q]Xi = 0} “the plane x\ — 0” simply. 

Figure 2 shows graphs of maximum values of |v| versus time t (left) and the distance from 
the maximum point of \v\ to the z -axis versus time t (right) for the swirl case, where four colors 
are used for the graphs of Re = 50,000 (red), 10,000 (green), 5,000 (blue) and 1,000 (purple). 
Figure 3 expresses corresponding graphs for the no swirl case. From figures 2 and 3 we can see 
that different phenomena appear in the swirl and no swirl cases. Typical phenomenon is the 
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Figure 1. Initial data for a — 1/8, left: swirl, right: no swirl 




Figure 2 . Graphs of maximum values of \v\ versus t (left) and the distance from 
the maximum point of \v\ to the 2 -axis versus t (right) for the swirl case with 
Reynolds numbers 50, 000 (red), 10, 000 (green), 5, 000 (blue) and 1, 000 (purple). 




Figure 3. Graphs of maximum values of \v\ versus t (left) and the distance from 
the maximum point of \v\ to the z-axis versus t (right) in the no swirl case with 
Reynolds numbers 50, 000 (red), 10, 000 (green), 5, 000 (blue) and 1, 000 (purple). 


drastic changes of the distance from the maximum point to the z-axis around t = 0.35 and 1.0 
in the swirl case, which does not appear in the no swirl case. The difference becomes clearer in a 
higher Reynolds number flow. In order to understand the difference we display time evolutions 
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Figure 4. Time evolution of the velocity magnitude |^| on the plane x\ — 0 
in the no swirl case with Re — 50,000. t — 0.1 (top left), 0.3 (top right), 0.7 
(middle left), 1.1 (middle right), 1.7 (bottom left) and 3.0 (bottom right). 

of velocity on the plane x\ — 0 for both cases with Re — 50, 000 in figures 4, 5 and 6. Figures 4 
and 5 show the time evolutions of \v\ on the plane for the no swirl and swirl cases, respectively, 
and figure 6 exhibits the time evolution of the axial velocity u z (^-component of v) in the swirl 
case. They imply that the flow dissipates straightforwardly as t increases in the no swirl case 
and that an interesting flow structure appears near the z-axis in the swirl case. 

From figures 2, 5 and 6 it is observed that a downward flow arises near the z-axis around 
t — 0.3, that the maximum value of |t| is attained near the z-axis and the lower boundary 
around the same time (t — 0.3), that a new upward flow arises near the z-axis around t = 1.3, 
and that the velocity attains its maximum value near the z-axis and the upper boundary around 
the same time (t = 1.3). 


Remark 2.4. The scale may be different in each figures. Red region represents high magnitude. 
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Velocity Magnitude 
6.185e-01 1 1.2 1.4 1.732e+00 


Velocity Magnitude 
6.185e-01 1 1.2 1.4 1.732e+00 

^^1 I 



Velocity Magnitude 
6.185e-01 1 1.2 1.4 1.7320+00 

I i I 


Velocity Magnitude 


6.185©-01 1 1.2 1.4 1.7320+00 



Figure 5. Time evolution of the velocity magnitude |p| on the plane x\ — 0 in 
the swirl case with Re = 50,000. t — 0.1 (top left), 0.3 (top right), 1.0 (middle 
left), 1.3 (middle right), 2.3 (bottom left) and 3.0 (bottom right). 


3. A MATHEMATICAL SUPPORT OF “INCREASING VELOCITY NEAR THE AXIS OF SYMMETRY 

AND THE LOWER BOUNDARY” 

When someone observes unstabilizing effects (such as increasing velocity in the time variation), 
it had better to compare with the stabilizing effects (such as regularity results mensioned in 
the introduction). In the numerical computation in the previous section, Giga-Hsu-Maekawa’s 
criterion (G-H-M’s criterion in the following) was in mind. G-H-M’s criterion (see A3.If) also 
( Giga et al . I 201 I Theorem 1.3.)) says that there are no type I blow-up solutions to (II.ip 


under a continuous alignment condition on the vorticity direction in the half space with no-slip 
boundary condition. More precisely, if the vorticity direction is controlled in some sense at the 
places where the absolute value of the vorticity is large, the flow is regular. Conversely, if we 
can construct solutions which do not satisfy such kind of continuous conditions at those places, 
we have more chances to find out the clue for possible blow-up solutions. Even in this paper 
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Figure 6 . Time evolution of the axial velocity u z on the plane x\ — 0 in the 
swirl case with Re — 50,000. t = 0.1 (top left), 0.3 (top right), 1.0 (middle left), 
1.3 (middle right), 2.3 (bottom left) and 3.0 (bottom right). Note that red and 
blue colors represent the positive and negative values in this figure. 


we could not construct a blow up solution, this kind of thought might help us to explain or to 
predict some phenomenons with unstabilizing effects (such as increasing velocity phenomenons 
in this paper). Note that in our numerical computation, even we use the cylindrical domain 
instead of the half space, the behavior near the center of the lower (and upper) boundary is 
expected to be similar to the boundary of the half space under the same boundary condition. 
Now we give a mathematical support why the maximum value is increasing near the saddle point 
on the boundary. Let u = uqcq + u r e r and u 1 - = u r eo — u$e r . We can calculate the vorticity 
on the boundary as uo — d z u ■ L . Note that d z u r (r, 0) -A 0 as r -A 0 due to the symmetry. We 
can show that the direction of the vorticity uj/\uj\ is not continuous at the saddle point. In this 
case we need to assume B z uq ^ 0 (on the boundary) and d r UQ ^ 0 (on the axis) near the saddle 
point. These nonzero conditions may express “shear flow effect by the swirl”. The vorticity 
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* Velocity Magnitude 

6.185e-01 1 1.2 1.4 1.732e+00 



Vorticity_Omega Magnitude 
1.580e-02 20 30 405.444e+01 



Vorticity_Omega Magnitude 
1.580e-02 20 30 405.444e+01 

am m"wmwmimM 


Figure 7. Contours on the plane x\ — 0 of the velocity magnitude |p| (first for 
t — 0.4, second for t = 1.4) and the vorticity magnitude |o;| (third for t — 0.4, 
fourth for t — 1.4) in the swirl case with Re — 50, 000. 
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Figure 8. Graphs of the quantities \uj\ (top) and the three components of £ 
(bottom) at t = 0.4 on the line parallel to the z -axis through the point 
(0,0.05,-0.125). 



Figure 9. Graphs of the three components of £ at t — 0.4 on the line parallel 
to the a^-axis through the point (0,0.05, —0.125). 


along the axis is expressed as uo — 2 d r UQe z . Thus uj/\uj\ along the z-axis and on the boundary 
is not continuous. Although, rigorously, it does not break the following continuous alignment 
condition in G-H-M’s criterion: 


(3.1) 


UJ , . UJ . 

— (t,x) - r(t,y) 

\UJ\ \UJ\ 


<p(\x-y |), 


where p(\x — y\) is any modulus continuous function, since the vorticity is zero at the saddle 
point. Nevertheless, by the numerical compuation, we see that there are high vorticity region 
near the saddle point and the boundary. We observe \uj\ and three components of £ — uj/\uj\ 
under a line which contains a point (0, 0.05, —0.125) on the boundary and parallels to z-axis at 
t = 0.4, see figure 8 (and see figure 9 for the case that the line parallels to a^-axis). The bottom- 
axis in figure 8 represents the distance from the lower boundary while the bottom-axis in figure 
9 represents the distance from the saddle point. Note that the coordinate of the saddle point is 
(0,0,-0.125) and the axial velocity u z — 0 on the boundary. At the point (0,0.05,-0.125) on 
the boundary, magnitude of the vorticity \uj\ is 60, and it clearly attains maximum value along 
the line in figure 8. Moreover uj/\uj\ is also drastically changing (highly oscillating) near the 
boundary (see figure 8 and figure 9). It shows some kind of unstability of £ at the location of 
the large vorticity magnitude. 
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4. Conclusion 


From literatures, it is natural to consider that the hyperbolic flow with swirl and saddle 
point on the boundary might be the key structure of flow and probable place for unstability 
effects occur near the no-slip flat boundary. We showed the clear structure for the axisymmetric 
hyperbolic flow with swirl and observed the following phenomenons which are distinctly different 
from those without swirl: (1) The distance between the maximum point of the velocity and 
the z-axis is drastically changing around some time we called it turning point. (2) The velocity 
increases and obtains its extreme value (maximum) near the axis of symmetry and the boundary 
when time is close to the turning points. By comparing these with the studies of tornadoes, it 
might help our understanding about the behavior of the velocity of wind near the ground which 
is very significant in the researches of tornadoes for reducing the damage cause by tornadoes 
or similar phenomenons. (3) The downward flow near z-axis is observed. (See figure 6 the 
time variation diagram.) The downward wind inside the core of real tornado is also observed 
in two-celled vortex structure of the studies of numerical simulations for time-averaged velocity, 


see (llshihara et all 2011, figure 4(b)). By comparing our observation with studies of tornadoes 


might enhance our understanding about the behavior inside the core of tornado for high swirl 
ratio. Furthermore, we also observed the properties of vorticity and its directions (see Section 3) 
which seems few in the literature of the numerical simulations of tornadoes. Those phenomenons 
might be the clues for helping our understanding on the behavior near the saddle point and near 
the boundary. And our numerical approach might be helpful for other boundary shapes and 
other equations related to fluid mechanics in the future work. It might be useful for the studies 
of tornadoes arise or pass by a different landforms instead of flat plane. 
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Appendix A. 

In this section, we present additional numerical results in order to show two things. The one 
is that the dependency of the numerical results in Section 2 on the discretization parameters h 
and r is qualitatively small, and the other is the reason of the choice of the range of z in 12, i.e., 
—a < z < 4 a. 

Firstly, we consider the former thing. We compute the no swirl and swirl cases for Re — 
50, 000, 10, 000 and 5, 000 by the stabilized Lagrange-Galerkin scheme with a coarse mesh and a 
large time increment. The maximum, minimum and average mesh sizes of the coarse mesh are 
2.48 x 1CT 2 , 1.93 x 10 3 and 1.20 x 10 2 , respectively, where the strategy of the mesh generation 
is the same, i.e., the mesh size around z- axis is smaller than that of other part. The time 
increment is set as r = 1.66 x 10 -2 . In the following we call the numerical results in Section 2 
“results A” and the numerical results by the coarse mesh and the large time increment “results 
B”. We compare results B with results A for Re = 50,000, 10,000 and 5,000. Figure 10 shows 
graphs of maximum values of |t| versus t (left) and the distance from the maximum point of |t| 
to the 2 ?-axis versus t (right) in the swirl case for Re — 50, 000 (top), 10, 000 (middle) and 5, 000 
(bottom), where red and green colors are employed for results A and results B, respectively. 
Figure 11 displays corresponding graphs in the no swirl case for Re = 50,000. From figures 10 
and 11, we can see that the two graphs in each figure are qualitatively similar, while there is 
difference quantitatively. Here, the graphs for Re = 10, 000 and 5, 000 in the no swirl case are 
omitted, since they are also similar qualitatively. We display additional information of results 
B in figure 12, which shows cross-sections of |t| (top) and |cj| (bottom) at t — 0.4 (left) and 
1.4 (right), for Re = 50,000 in the swirl case. By comparing figure 7 (results A) with figure 
12 (results B) we can see that the behaviors are almost the same, although the magnitudes are 
slightly different. These results imply the former thing, i.e., the dependency of the numerical 
results in Section 2 on the discretization parameters h and r is qualitatively small. 

Now we consider the latter thing. Let Q! be a domain defined by Q' := {{x\,X 2 ,z) E M 3 : 
—5a/2 < z < 5a/2, \/x\ + x\ < 1}. We note that centers of ft and are (0, 0, 3a/2) and 
the origin, respectively. We compute the swirl case in the domain IF for Re = 50, 000 by the 
stabilized Lagrange-Galerkin scheme with a coarse mesh and a large time increment, where the 
initial data is set by the same functions in (|2.3|) - (l2.6p . The mesh sizes and the time increment 
are the same as used in results B. (The mesh is generated by a translation of the mesh for results 
B.) We call the numerical results “results C”. We compare results C with results B in figure 13, 
which shows graphs of maximum values of |t| versus t (left) and the distance from the maximum 
point of \v\ to the z- axis versus t (right) in the swirl case for Re = 50, 000, where red and green 
colors are employed for results B and results C, respectively. It is observed that the “increasing 
velocity phenomenon” in results B is clearer than that in results C, and that drastic changes of 
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Figure 10. Graphs of maximum values of |p| versus t (left) and the distance 
from the maximum point of |i;| to the 2 -axis versus t (right) in the swirl case for 
Re = 50,000 (top), 10,000 (middle) and 5,000 (bottom), where red and green 
colors are employed for results A and results B, respectively. 

the distance from the maximum point of |p| to the 2 -axis are observed in both results. Hence, ft 
is better than ft' in order to see both of the “increasing velocity phenomenon” and the “drastic 
changes”. These results explain the latter thing, i.e., the reason of the choice of the range of z 
in ft, —a < z < 4 a. 


References 

Bourgain, J. & Li, D. 2013 Strong ill-posedness of the incompressible Euler equations in 
borderline Sobolev spaces. To appear in Inventione Math, preprint arXiv:1307.7090 [math.AP] 


































A LOCAL ANALYSIS OF THE AXI-SYMMETRIC NAVIER-STOKES FLOW 


15 


co 
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for Re = 50, 000, where red and green colors are employed for results A and 
results B, respectively. 
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